arXiv: 1505.01084v2 [math.PR] 17 Jul 2015 


CENTRAL LIMIT THEOREM UNDER UNCERTAIN LINEAR 

TRANSFORMATIONS 


DMITRY B. ROKHLIN 

Abstract. We prove a variant of the central limit theorem (CLT) for a sequence of 
i.i.d. random variables £j, perturbed by a stochastic sequence of linear transforma¬ 
tions Aj, representing the model uncertainty. The limit, corresponding to a ’’worst” 
sequence Aj , is expressed in terms of the viscosity solution of the G-heat equation. In 
the context of the CLT under sublinear expectations this nonlinear parabolic equation 
appeared previously in the papers of S. Peng. Our proof is based on the technique 
of half-relaxed limits from the theory of approximation schemes for fully nonlinear 
partial differential equations. 


1. Problem formulation 

Consider a sequence of i.i.d. d- dimensional random variables Q = (£j)f =1 . 

Denote by £ a random variable distributed as £j, and assume that 

Ef = 0, (E(f ?'))', =1 = /, (1.1) 

where I is the identity matrix. By the classical central limit theorem (CLT), for any 
bounded continuous function / : t —y K. and a fixed d x d matrix A, we have 

lim E f (it, r ) = 

where 77 has the standard d- dimensional normal law. Note, that for given / the limit 
depends only on the covariance matrix AA T of 

In this paper we consider the case where A is not known exactly, and can change 
dynamically within a prescribed set. This is a simple example of a probability model 
under uncertainty. The extension of the CLT, obtained below, looks similar to Peng’s 
CLT under sublinear expectations: mm- However, our problem formulation, as well 
as the proof, do not involve the nonlinear expectations theory in any way. On the other 
hand, similarly to Peng’s approach, the key role is played by the viscosity solutions 
theory. 

Consider a filtered probability space ( 0 , J^", P, (J^)°h 0 ) and an adapted sequence 
1 °f d-dimensional random variables such that is independent from \ and 
satisfy (11.11) . Denote by (resp., E>d) the set of d x d matrices (resp., symmetric 
matrices). Let (A,)“L 0 be an adapted sequence with values in a compact set A e M^. 
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The process (Aj, &j)JL Q , where is j > 1 may be wider than ct(£l,..., £j), is chosen 
by the “nature”, and represents the “Knightian uncertainty”. 

Our goal is to describe the limit 


Jgf : = 


lim 


sup 


E / 



( 1 . 2 ) 


where / is a bounded continuous function on R d , A™~ 1 = (Aj)™~l and 'QGf 1 = 
{A™~ 1 \Aj : 0 n- A and Aj is ^-measurable}. Besides the theoretical interest, such 
quantities are useful, e.g., for measuring risk and option pricing under uncertain volatil¬ 
ity: see [I]. 

Our result (Theorem [1]) and its discussion are presented in the next section. The 
proof is deferred to Section [3] In Section [4] we discuss the relationship of our problem 
to the nonlinear expectations theory. 


2. Central limit theorem under uncertain linear transformations 


First we formulate our result and then give some comments. 


Theorem 1 . Let v : [0,1] xl d i-)R be the unique continuous viscosity solution of the 
nonlinear parabolic equation 

- v t (t, x) - - sup Tr (. AA T v xx (t , x)) = 0, (t, x ) £ [0,1) x R d , (2.1) 

2 AeA 


satisfying the terminal condition 

u(l, x) = f(x), x £ R d . 


( 2 . 2 ) 


Then 2zf = n(0, 0). 


The equation (12. ip can be written in the form 

—v t — G(y xx ) = 0, where G(S) = -supTr ( AA T S ) , S £ § d . 

2 ag a 

By v xx = [v x r x i)f l=1 we denote the Hessian matrix. 

Let us recall the definitions of viscosity semisolutions: see, e.g., [8j. Put Q° = 
[0,1) x R d , Q = [0,1] x R d and denote by Cf(M d+1 ) the set of test functions, whose 
derivatives up to the second order, are continuous and bounded. A bounded upper 
semicontinous (use) (resp., lower semicontinuous (lsc)) function u : Q t —> M, is called a 
viscosity subsolution (resp., supersolution) of the problem (12.11) . (j2.2j) if u(l,x) < f(x) 
(resp., u(l,x) > f(x)) on R d and for any ip £ C%( M d+1 ), (t,x) £ Q° such that (t,x) is 
the strict global maximum (resp., minimum) point of u — ip on Q°, the inequality 

-<p t (i,x) - G (p xx (t, x)) < 0 (resp., > 0) 

holds true. As is well known, in these definitions one can equivalently assume that the 
extremum is local or not strict. 

A bounded continuous function u : Q H» R is called a viscosity solution of (12.11) , 
(12.2j) . if it is a viscosity sub- and supersolution. 
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Remark 1. The uniqueness of a bounded continuous viscosity solution of (12.11) . (12.21) 
is well known: see, e.g., [TU] (Theorem 3.5). However, to prove Theorem |T] we need 
a more subtle result. Note, that we require the equation (12. ip to be satisfied in the 
viscosity sense at the lower boundary of Q. So, by the accessibility theorem of [6], for 
an use viscosity subsolution u and a lsc viscosity supersolntion w of (12.ip we have 


u( 0 ,t)= lirnsup u(t,y ); iu(0,t) 

(t,y)€(0,l)xR ti , 

t—» 0 , 


liminf w(t,y). 

(f,?/)e(o,i)xR d , 
t—>0,y—tx 


Using this fact we may apply the comparison result of [2j (Theorem 1) and conclude 
that u < w on Q. 


Remark 2. Consider the system of stochastic differential equations 

dY s = A s dW s , s G [t, 1], Y t = x, (2.3) 

where W is a standard d -dimensional Brownian motion. Denote by At the set of 
stochastic process (H s ) se [ ti ij, which are progressively measurable with respect to the 
minimal augmented filtration, generated by W (see Chap. 1 of [5]), and take values in 
A. From the theory of stochastic optimal control (see mm) we know that the value 
function 

V{t,x)= sup E f(Y 1 ) (2.4) 

AeAt 

is a continuous viscosity solution of the Hamilton-Jacobi-Bellman equation (12.ip . Tak¬ 
ing (t, x) = (0,0), from Theorem [T] we obtain the stochastic control representation of 
the limit (II , 2 p : 

^ = U(o, 0 ) = sup E f([ dW 8 ) . 

AeA 0 \J o / 

Remark 3. The Euler scheme 

b j/n+l/n b j/ n T Aj/ n (fUj/n+l/n j 0, . . . , TL 1 

for ( 12 .3p gives the approximation 

(n-l A \ 

m0)« sup Ef V j/n ja +1 , r h+l = V^(W j/n+l/n - W j/n ) 

A j/n eA \^ 0 Vn ) 

for the value function (12.4p . similar to the expression in (1 1.2 f) . However, in contrast 
to £j, are normal. From Theorem |T] it follows that to get a correct approximation of 
(12.4p . one can take in the Euler scheme instead of ry any i.i.d. random vectors 
satisfying ( 11 . 11 ) . 

Remark 4. Equation (12. ip is called G-heat equation: [14], [13). ft was used by S.Peng 
for the description of the G-normal distribution. Moreover, the representation of the 
same form, as given in Theorem [lj appeared in the CLT under sublinear expectations: 
see [421 QH . 

Remark 5. Similar to the classical case, for fixed /, the limit AY depends only on the 
set {AA T : A G A} of possible covariance matrices of A £. We stress that SY does not 
depend on the choice of filtration D c(^i,..., £j). 
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To prove Theorem [H let us introduce the state variables Xj by 



X j+1 = Xj + ^±, X s =x, 
y'n 

j = s,...,n- 1 . 

(2.5) 

Denote the solution of (]2.5[) by X S,X,A and consider the value functions 



Vn0-,X) = f{x), 


( 2 . 6 ) 


V n (s/n,x) = sup E/ (X s ^ xA ) , 
as-Rqi ™- 1 

s — 0 ,..., n — 1 . 

(2.7) 

Clearly, 

& = lim sup E/ (X°’°’ A ) 

n —*°° Tq — 1 — 1 

= lim v n (0, 0 ). 

n— >• oo 

( 2 . 8 ) 


Our goal is to prove that u n (0, 0) —> u(0,0), where v defined in Theorem [TJ We 
apply the half-relaxed limits technique of [3], which became standard in the theory of 
approximation schemes for fully nonlinear second order elliptic and parabolic equations 
after the seminal paper [4]. Namely, we construct the half-relaxed limits v < v of 
v n , and prove that they are viscosity semisolutions of (12. ip . (12.2p . Then we use the 
comparison result, mentioned in Remark [TJ to prove the opposite inequality: v>v. It 
follows that v = v = v and v n —> v. 

Note, that the strategy of HUE] is different and based on the following fact: 

lim sup Eu(1,X°’ M ) =u(0,0). (2.9) 

The direct proof of (12. 9 p requires the interior Holder regularity of v, which is guaranteed 
in the strong parabolic case, and Lipshitz continuity of /. The general case requires 
perturbations and approximations: see [14] . 

3. The proof of Theorem [T] 

The discrete time stochastic control problem (12.5p . (12.6j) . (12.71) looks quite standard. 
However, to model Knightian uncertainty, we consider the class A n of “open-loop” 
strategies Aj, adapted to an arbitrary filtration D cr(£i,... ,£j) (cf. P2])- So, we 
prefer to give a direct proof of the dynamic programming principle, instead of trying 
to find an appropriate reference. 

Lemma 1. Put t™ = j/n,j = 0,... ,n. The value functions v n are continuous in the 
state variable x and uniformly bounded: 

\v n {t],x)\ < sup \f(x)\. (3.1) 

rrGR d 

Moreover, they satisfy the recurrence relations 

v n ( l,x) = f(x), 

v n (tf, x) = sup E v n {t n j+l , x + A£/y/n), j = 0 ,..., n - 1 . 

Te A 


(3.2) 

(3.3) 
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Proof. Consider the sequence v n , defined by (I3.2j) . (13.3ft . Clearly, u n (l, •) is continuous 
and bounded. Assume that v n (t™ +1 , ■) has the same properties. Then Ff(x, A) = 
Ev n (t" +V x + Af/y/n) is continuous in (x,A) by the dominated convergence theorem, 
and the function 

v n (t ”, x ) = sup Ff(x, A) 
aga 

is bounded and continuous in x by the compactness of A. 

Furthermore, by the Dubins-Savage type measurable selection result [18j , [20] (The¬ 
orem 5.3.1), there exists a Borel measurable function : R d H> A such that 

v n (tf,x) = Ff(x, $j(x)) = E v n (t] +1 ,x + $j(x)£/y/n), x G R d . (3.4) 

Define the process X by the recurrence relations 


Xj + i — Xj + $j(Xj)£j + i/ \fn, X s — x, 


(3.5) 


and put A* = &j(Xj). The process (13.5p can be written as X S,X,A *. From (13.41) we get 


Vn(t],x;-’’ A ')= I V n (t] +1 ,X’^ + <S>j(X 


rS,x,A* 


s,x,A* 


)zj\fn) dPc(dz), 


where is the distribution of £. Taking the expectation and using the independence 
of X*' x ' A and ^j+i, we obtain 


E v n {i 


n \r 
j ? 


s,x,A* 


= E v, 


.(<”+ i,xp A ' + tAxp A ')(ni/V^) = E«„((” +1 ,x*;y). 


Thus, 

v n K,x) = Ev n K,Xr- A ‘) = Ev n (CX’^ A ') = E f(X‘-*- A '). (3.6) 

On the other hand, v n (tf,x) > E v n (tf +1 ,x + A^/^/n), AeA, and for any sequence 
(Aj)jfs e 2t” _1 we have 


A'"" 4 ) > 


VnXj+i.Xj A ' A + Ajzjy/n) dpdd:). 


Taking the expectation, and using the independence of (Xj' x ' A , Aj) and £j+i, we obtain 


E v n (tf,Xp A ) > Ev n {tf +1 , Xj ,x,A + Aj£ j+1 /y/n) 


s,x,A 


E v n (t 


n 

j+ 1’ 



Hence, 


VnK.x) = Ev n (t:,Xr A ) > Ev„(C,X"' A ) = Ef(X A ’’ A ) 


(3.7) 


for any {Aj) n -Z] G Combining (j3.6[) and (13.71) . we conclude that the function, 

defined by the recurrence relations (13.21) . (13.31) . is the same as the function (j2.6j) . (12.71) . 
The inequality (13.11) follows from (12.71) . □ 


As follows from the proof, an “optimal strategy” A* of the nature uses only the 
information on the current state Xj, although the available information may be 
much richer (see [19] for a similar conclusion). 

Consider a closed set U G M m and a sequence U n C U of its closed subsets such that 
for any x G U there exist exists a sequence Xk G U nk , G N converging to x. For 
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a uniformly bounded sequence of continuous functions u n : U n (->• M, \u n \ < M define 
Barles-Perthame type half-relaxed limits u,u:U KRas follows 

u(x) = inf{limM nfc (x fc ) : x k G U nk , x k ->■ x and u nk (x k ) converges}; 
u(x) = sup{limu nfe (x fc ) : x k G U nk , x k ->• x and u nk {x k ) converges}. 

It follows from the definitions that there exist sequences n k G N, x k G U Uk such that 
x k ->■ x, u nk (x k ) ->• w(x) (resp., u nk (x k ) -> u(x)). 

The proofs of the next two lemmas follow the argumentation of [2| (Chap. V, Lemmas 
1.5, 1.6). 

Lemma 2. u is Isc, u is use. 

Proof. If u is not lsc at x G I/, then there exist 5 > 0, J G N and a sequence y 3 G U, 
Dj —» x such that w(r/j) < u(x) — 5 for j > J. By the definition of u for each j there 
exist rij gN, x 3 G U nj , \x 3 — y 3 \ < 1 /j such that 

u nj {xj) < u(yj) + 5/2. 

Thus, u nj (xj ) < u(x) —5/2, j > J in contradiction with the definition of u. 

The case of u is considered in the same way. □ 

Lemma 3. Let tp G x G U. Ifxis the strict global minimum (resp., maxi¬ 

mum) point of u — ip (resp., u — ip) on U, then there exist sequences n k G N, y k G U n . k 
such that y k —* x, u nk (y k ) —> u(x) (resp., u nk (y k ) —> u(x)), and a test function if such 
that ip x (x) = <p x {x), ifxx(x) = (Pxx{x) and y k is a global minimum (resp., maximum) 
point of u nk - if on U rik . 

Proof. We consider the case of u. Let a sequence x k — > x, x k G U nk be such that 
u n k (x nk ) —> u(x), and let y k be a minimum point of u nk — p on the set U nk D B i, 
B\ = {y '■ \y — x\ < 1}. Then 

(u nk - p)(y k ) < ( Un k - <p)(x k ) (3.8) 

for sufficiently large k. Passing, if necessary, to a subsequence, we may assume that 
y k —* y G U, u nk (y k ) —> z. Thus, by (13.81) and the definition of u, we get 

u(y) - <p(y) <Z- pi/y) < u(x) - p(x). 

It follows that y = x, z = u(x), since x is the strict minimum point of u — p on U. 
Furthermore, y k —* x is a local minimum point of u nk — p on U nk , since y k lies inside 
the ball Bi for sufficiently large k. 

Let x G C/ be a function such that y(x) = 0, |x — x| < 1/2, x(x) = 1, |x — x| > 1. 
Then there exist M' > 0 such that y k is a global minimum point of u n . k (x ) — p{x) + 
M'x(x) on U nk for k large enough. The sequence y nk and the test function if = p — M'x 
have the desired properties. □ 

In the theory of approximation schemes the following relation is known as a consis¬ 
tency condition: see [lij. 

Lemma 4. Let G C/(M x M d ) and (t n ,x n ) —> (t,x). Then 

lirn n sup E (p(t n + 1/ n, x n + A£/^/n) - p(t n ,x n )) = p t (t,x ) + G (p xx (t,x)) . 

n^oo AgA 
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Proof. By Taylor’s formula we get 

<p{tn + l/n, x n + At,/y/n) - <p{t n ,x n ) = <p(t n + 1 /n,x n + At,/y/n) 

-<p{t n , x n + Af/y/n) + <p(t n , x n + At/y/n) - ip(t n , x n ) 

1 i i 

Ptifni Xn T A // y/n) T ~^Px(fni Xn)A/ T — ( Pxxifni X n )A^ • Af, (3.9) 

n y/n in 

where t n — t n + a n /n, x n = x n + f3 n Af/y/n, a n , /3 n G [0,1]. Note, that the sum of the 
first and third terms in the last line of (13. 9 p is a (£)-measurable, and the expectation 
of the second term is 0 since E£ = 0. 

Using (13.91) . we get the inequality 


n sup E (<p(t n + l/n,x n + A£/y/n) - ip(t n ,x n )) - (<p t (t,x) + G (ip xx (t,x))) 

Ae A 


< sup 
Ae A 


^-Ptifni X n T A/ / y/n) T ^ E ( iPxxifni x/)A(f • 


-(Pt(t,x) - -Tr [AA T ip xx {t,x)) 


which yields the result by the dominated convergence theorem. 


□ 


To apply Lemmas El [3] in our case, put m — d + 1, U — Q — [0,1] x M d , U n = Q n = 
U;=o {j/ n } x ® d ■> and denote by v, v the half-relaxed limits of v n . 

Let (t, x) G Q be the strict global minimum point of v — tp on Q for a test function 
ip G Take sequences y k = (tj( k ),x k ) G Q nk —> (t,x), v nk and a function if, 

given by Lemma [3] If t < 1, we may assume that < 1, that is, j(k) < n. By the 
recurrence relation (13.31) we have 

v nk (tj(k),x k ) = sup E v nk {tj( k \ + l/n k ,x k + A£/y/n k ). 

Ae A 

The inequality (v nk - if)(t,x) > (v nk - if)(t j{k) ,x k ), (t,x) G Q nk implies that 

0 > sup E if[t j(k) + l/n k , x k + A£/y/n k ) - ip(tj(k),x k ), 

Ae A 

and Lemma E] gives the inequality 

0 < -cp t (t,x) - G (<p xx (t,x)) , (3.10) 

since the derivatives of ip and if, up to the second order, coincide at (t,x). 

Now assume that t = 1. Clearly, u(1,t) < f(x). If n(l,T) < f(x) then tj( k ) < 1 for 
sufficiently large k (since v nk (l,x k ) = f(x k ) converges to f(x)). So, we again obtain 
the inequality (13. lOj) as above. However, this is impossible. Indeed, we can change ip 
to (p = ip — c( 1 — t), c > 0 in this inequality, since ( 1 ,T) is still the global minimum 
point of v — (p\ 

0 < -c - ipt - G (ip xx (t, x)) , for any c > 0. 

This contradiction shows that vfl,x) = f(x), x G M d . 

Thus, we have proved that v is a viscosity supersolution of (12. ID . (12.21) . In the same 
way one can prove that v is a viscosity subsolution of (12 .1 p . (12.2h . By the comparison 
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result, mentioned in Remark [U we have v > v. The opposite inequality follows from 
the definition of v, v. Therefore, v — v — v is the unique continuous viscosity solution 

of (im m- 

Finally, form the definition of v,v we see that 

u(0, 0) = u(0, 0) < liminf u„(0, 0) < lim sup u n (0, 0) < u(0, 0) = u(0, 0). 

n—>oo n—>-oo 

In view of (12.8p . this hnishes the proof of Theorem 1. 


4. On the relationship with the sublinear expectations framework 

Recall (see, e.g., m that a sublinear expectation space is a triple (0,%,E), where 
Dis a set, V. is a linear space of real valued functions defined on D, and E is a sublinear 
functional on "H. It is assumed that "H contains constants and |A"| G W for X G H. 
Moreover, "H should be invariant with respect to some functional transformations. The 
most standard assumption is the following: 

<p( Xi, ...,X n )eH for x lt ..., X n eu, ipe C lMp (W l ), 

where C'i i Lj P (R ri ) is the linear space of functions <p, satisfying the inequalities 
\<p(x) - (p(y)\ < C(1 + \x\ m + \y\ m )\x - yl x,y el n , 
with C , m, depending on ip. 

Sublinear expectation EAaI satisfies the following conditions: 

(i) EX < E Y, X < Y, 

(ii) Ec = c for a constant c G R, 

(iii) E(X + Y) < EX + E Y, 

(iv) E(AX) = AEX, A > 0. 

To put our problem in this context, assume that the random variables £* have finite 
moments of all orders. Consider the space of sequences = {(xj )“ 1 : x^ G M d }, and 
introduce the space of random variables % as follows: % = U '^ > =1 'H n , where T~L n is the 
space of continuous functions X = ijj(xi,...,x n ) of polynomial growth. Define the 
sublinear expectation by the formula 

ELY = SUp EV>(A)£l, ■ ■ ■ , 

It is easy to see that the triple (D, %. E), dehned in this way, is a sublinear expectation 
space. 

Denote by X t the projection mappings: X^x) = x % . For p G Ci } L ip (R dxn ) put 

Tp{xi, .. .,x n - 1 ) = E <p(x !,... ,x n _i,X n ) = sup Ep{xh ... ,x„_i ,A^ n ). 

Ae A 
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Using the result of [16] (Theorem 14.60) on the interchange of maximization and ex¬ 
pectation operations, it is not difficult to show that 

E<p(Xi,...,X n ) = sup E<p(A 0 fi,... ,A n _^ n ) 

= sup E^(A 0 £i, • • •, A n _ 2 in~i) = E^(Xi,..., X n _i). 

A %~ 2 G %~ 2 

This means that X n is independent from (Xi, ..., X n _i) in the sense of sublincar 
expectations theory (see [I5j, Definition 3.10). 

Denote by m n ( A) the set of .^-measurable functions with values in A. The random 
variables X t have no mean uncertainty: 

E(±X;) = sup E(±Aj_i£j) = 0. 


Furthermore, for S G we obtain 

V.(SXi-Xi)= sup E(SAi-i£i ■ Ai-iti) = sup ETr (Aj^SA^) 

Ai_iGmi_i(A) Ai_iGmi_i(A) 

= sup Tr (A t SA). 
aga 


By Peng’s central limit theorem (see |I4) Theorem 5.1], [151 Theorem 3.3]) the sequence 
n -i /2 Xi converges in law to a G-normal random vector Y: 


E / 



'n—l 


sup 


U £ 


o=o 


Aj£j+ i 
\fn 


-> i/on, 


/ e C b (R) 


with 

G{S) = ^E (SXi ■ Xi ) = i sup Tr (A T SA). 

4 2 agA 

By the definition of the G-normal distribution, we have E f(Y) = u(0, 0), where u is 
the viscosity solution of the G-heat equation (12.ip with the boundary condition (12.21 . 
Hence, for £, with hnite moments of all orders, Theorem 1 follows from Peng’s CLT. 

As is mentioned in [15] (Remark 3.8) the same result can be proved under the as¬ 
sumption E|JXi| 2+<5 < oo for some 5 > 0, instead of E|A 7 i | n < oo for all n (it is not known 
whether one can take 5 = 0). In our case it corresponds to the condition E|^| 2+<5 < oo, 
which is still superfluous. The reason for the appearance of this assumption we see in 
the lack of dominated convergence theorem in a general sublincar expectation space. 
We also refer to ra. where Peng’s approach was applied to the one-dimensional prob¬ 
lem with variance uncertainty. It was shown that, written in the classical terms, this 
approach allows to prove the CLT without unnatural assumptions, even in the case of 
non-identically distributed independent random variables. 

Conversely, one can try to prove Peng’s CLT by the methods of the present paper. 
Let (W)^! be i.i.d. random variables under a sublinear expectation E, and let / be 
a bounded Lipshitz continuous function. The definition of independence, applied to 
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bounded Lipshitz continuous functions 


Vn{tj, x)=Eflx + Xi ) > *]= i/ n ’ 


*=3 


gives the recurrence relation 
v n (t],x) = I 


Xj 1 


E / D + fffe + TR E Y 


n 1 


i=j +1 


= Ev n (tj +1 , x + Xj/n 1/2 ) 


corresponding to the dynamic programming principle, considered in Lemma [T] How¬ 
ever, to follow subsequent reasoning, one needs additional assumptions, like monotone 
continuity (or Fatou) property: see [7j. Indeed, we have used the dominated conver¬ 
gence theorem which is not true in a general sublinear expectation space. 


References 

[1] M. Avellaneda, A. Levy, and A. Paras. Pricing and hedging derivative securities in markets with 
uncertain volatilities. Appl. Math. Finance, 2(2):73-88, 1995. 

[2] M. Bardi and I. Capuzzo-Dolcetta. Optimal control and viscosity solutions of Hamilton- Jacobi- 
Bellman equations. Birkhauser, Boston, 1997. 

[3] G. Barles and B. Perthame. Exit time problems in optimal control and vanishing viscosity method. 
SIAM J. Control Optim., 26(5): 1133-1148, 1988. 

[4] G. Barles and P. E. Souganidis. Convergence of approximation schemes for fully nonlinear second 
order equations. Asymptot. Anal, 4:271-283, 1991. 

[5] R.F. Bass. Stochastic processes. Cambridge University Press, Cambridge, 2011. 

[6] Y.-G. Chen, Y. Giga, and S. Goto. Remarks on viscosity solutions for evolution equations. Proc. 
Japan Acad., Ser. A, 67(10):323-328, 1991. 

[7] S. Cohen, S. Ji, and S. Peng. Sublinear expectations and martingales in discrete time. Preprint 
arXiv:1104.5390 [math.PR], 26 pages, 2011. 

[8] M. Crandall, H. Ishii, and P.-L. Lions. User’s guide to viscosity solutions of second-order partial 
differential equations. Bull. Amer. Math. Soc., 27:1-67, 1992. 

[9] J. Diehl, P.K. Friz, and H. Oberhauser. Regularity theory for rough partial differential equations 
and parabolic comparison revisited. In Stochastic Analysis and Applications 2014, pages 203-238. 
Springer, 2014. 

[10] M. Nisio. Stochastic control theory: dynamic programming principle. Probability Theory and 
Stochastic Modelling, 72. Springer, Tokyo, 2015. 

[11] S. Peng. G-expectation, G-Brownian motion and related stochastic calculus of Ito type. In Sto¬ 
chastic analysis and applications: The Abel Symposium 2005, pages 541 567. Springer, 2007. 

[12] S. Peng. Law of large numbers and central limit theorem under nonlinear expectations. Preprint 
arXiv:math/0702358 [math.PR], 8 pages, 2007. 

[13] S. Peng. Multi-dimensional G-Brownian motion and related stochastic calculus under G- 
expectation. Stock. Proc. Appl, 118(12):2223-2253, 2008. 

[14] S. Peng. A new central limit theorem under sublinear expectations. Preprint arXiv:0803.2656 
[math.PR], 25 pages, 2008. 

[15] S. Peng. Nonlinear expectations and stochastic calculus under uncertainty. Preprint 
arXiv:1002.4546 [math.PR], 149 pages, 2010. 

[16] R.T. Rockafellar and R.J.-B. Wets. Variational analysis. Springer-Verlag, Berlin, 2009. 

[17] D.B. Rokhlin. Central limit theorem under variance uncertainty. Preprint arXiv:1506.01551 
[math.PR], 9 pages, 2015. 

[18] M. Schal. A selection theorem for optimization problems. Arch. Math., 25(l):219-224, 1974. 






CENTRAL LIMIT THEOREM UNDER UNCERTAIN LINEAR TRANSFORMATIONS 


11 


[19] M. Sirbu. A note on the strong formulation of stochastic control problems with model uncertainty. 
Electron. Comm. Probab., 19:1-10, 2014. 

[20] S.M. Srivastava. A course on Borel sets. Springer-Verlag, New York, 1998. 

[21] N. Touzi. Optimal stochastic control, stochastic target problems, and backward SDE. Fields Insti¬ 
tute Monographs, 29. Springer, New York, 2013. 

Institute of Mathematics, Mechanics and Computer Sciences, Southern Federal 
University, Mil’chakova str., 8a, 344090, Rostov-on-Don, Russia 
E-mail address: rokhlin@math.rsu.ru 



